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Abstract 

In acoustics, ultrasonics and in electromagnetic wave propagation, the 
crossed medium can be often modelled by a linear invariant filter (LIF) 
which acts on a wide-sense stationary process. Its complex gain follows 
the Beer-Lambert law i.e is in the form exp [— az] where z is the thickness 
of the medium and ct depends on the frequency and on the medium prop- 
erties. This paper addresses a generalization for electromagnetic waves 
when the beam polarization has to be taken into account. In this case, 
we have to study the evolution of both components of the electric field 
(assumed orthogonal to the trajectory). We assume that each component 
at 2 is a linear function of both components at 0. New results are obtained 
modelling each piece of medium by four LIF. They lead to a great choice 
of possibilities in the medium modelling. Particular cases can be deduced 
from works of R. C. Jones on deterministic monochromatic light. 

keywords: linear filtering, polarization, Beer-Lambert law, random 
processes. 

1 Introduction 

1.1 The Beer-Lambert law 

The Beer-Lambert law (B.L law) states that some positive quantity A (0) at the 
input of some medium varies following the equation 

A(z) = A(0)e- az (1) 

where z is the covered distance and a is a parameter defined by the medium 
[5J. The equality (1) comes from the approximation 

A (z + dz) - A (z) fa -aA (z) dz 

which postulates that the evolution of A (z) on a small thickness dz of the 
medium is proportional to A(z) and to dz, with a coefficient a > which is 
defined by the medium. Then the differential equation A' (z) — —aA (z) which 
leads to (1) . A more general method starts from the equality 
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A (0) A(z + z') = A (z) A (z') (2) 



1 



whatever z, z'. Equivalently the quotient A (z + z') /A (z) depends only on z\ 
and any piece of the medium of length z has the same behavior. If we assume 
that A(z) is a continuous function on R + , the only solution of (2) is (1) for 
some aeC. 

If we take A (0) = e 1 " *, (1) becomes (1Z [a] and X [a] stand for the real and 
imaginary parts of a) 



A (z) = exp 



iu>o ( t X [a] I — zlZ [a] 

LUQ 



(3) 



This means that the monochromatic wave e lu)at is delayed by ^X [a] and weak- 
ened by exp[— zlZ [a\] when crossing a thickness z in the medium. 

Now we place ourselves from a signal processing perspective. We assume that 
a piece of medium of any thickness z is equivalent to a LIF (Linear Invariant 
Filter) T z with complex gain F z (lu) and that any piece of thickness z + z' has the 
behavior of two filters in series T z and T z > |15j . [5]. This means that whatever 
the frequency Lu/2ir 

F z+Z , (u) = F z (u) F z , (w) . (4) 

Obviously we take Fq = 1. What preceeds implies that for each lu it exists a 
complex a (lu) such that 

F z (lu) = e-"*M 

By definition F z (lu) e lu}t is the output of the filter T z when the input is e lujt . 
For such an input the power P z at the output is 

P z = e - 2 ^K")] ( 5 ) 

(5) summarizes the Beer-Lambert law for wave propagation through a contin- 
uous medium, used in acoustics, ultrasonics and electromagnetics, a (lu) gives 
the attenuation of the wave (by its real part) and the celerity of the wave (by 
its imaginary part). The Kramers-Kronig relation links the real and imaginary 
parts of F z (w) which constitute a pair of Hilbert transforms [17], [11]. (5) is 
matched to monochromatic waves. More generally when the stationary process 
Z = {Z (t) ,t € R} is the input of F z , the output Z z = {Z z (t) , t G K} verifies 
(E[..] stands for mathematical expectation or ensemble mean) 



Z z (t)\ 2 = / e-^WsoMdu 



P*=E 



where so (w) is the power spectral density of Zo (so (w) = 5 (lu — luq) for a unit 
power monochromatic wave at luq). The fact that measurements are generally 
performed by non-monochromatic waves lead to gaps with the Beer-Lambert 
law in the form (6) [2]. Also, B.L law can be untrue when powers are too 
high (which act on the medium properties) [T] or when the beam expands [10) . 
However, B.L law has many applications in physics, chemistry and medecine 

0, 0, a, na. 

1.2 A counter-example 

The Beer-Lambert law applies for light amplitude or power after crossing contin- 
uous permanent media which can be viewed as a set of filters in series. However 
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light or radar wave is not only denned by an amplitude [3J, [13] ■ Electric field 
W- = (E*,EJj) is a two-dimensional vector with respect to axes Ox and Oy at 
distance z of origin O and with components 

E* = {El (t) ,t E R},E z y = {El (t),t e M} . 

It is assumed that the beam propagates in the neighbourood of the axis Oz and 
that the field is orthogonal to this axis. The beam is polarized in the direction 
9 at z when 

E% (t) =A(t) cose 
E* (t) = A (i) sin# 

for some (real or complex) process A = {A (t) , t € M} . The beam is unpolarized 
at z when the cross-spectrum s z y (uj) of components verifies 

4y H = o 

whatever the basis Oxy. This implies the equality of power spectra. The wave 
is partially polarized in other cases. 

The power P z at z for a stationary wave is defined by 



P 2 =E 



E z x (t)\ 2 +E \E* v (t)\* . (6) 
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We know that media act upon polarization and then can influence measure- 
ments, for instance in the case of antennas (generally matched to a particular 
polarization state). B.L law is not available for behavior of a given component 
of the electric field except particular cases. For instance, take the wave 

E°(t) =e W0 \ E°(t) = 

which propagates in a medium which rotates the beam by angle proportional to 
thickness z. We have (c is the celerity in the medium) 

El(t) = e l " a{t - z/c) cos [za(cj )\ 

If we measure E£ (t) the medium is equivalent to a filter of complex gain 

F z (w) = e ~ tulz/c{ui) cos [za (w)] 

when a particular direction is chosen (F z (uj) = e ~ tulz / c sin [za (uj)] for the or- 
thogonal direction and c can depend on uj). The term a (uj) takes into account 
a possible dependency of the rotation angle with the frequency. For an antenna 
which selects the component in a given direction at distance z, the B.L law is 
not verified. 



1.3 The two-dimensional case 

The B.L law is established for one-component monochromatic waves (i.e for 
waves defined by only one quantity depending on the time t and on the space 
coordinate z). They are time- functions in the form ae 1 " * (a G C depends on z). 
Quasi-monochromatic waves can be defined as random processes in the form 

E z (t) = A z (t) e tUot 
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where A z = {A z (t) ,t e R} is stationary with a baseband spectrum which can- 
cels outside (— ui,u!i) with ui/wo <C 1. The spectral band of E z is included in 
the interval (wo ~ ^1,^0 + ^1) ■ 

In the two-components case (i.e for waves defined by two quantities E z (t) , E z (t) 
depending on the time t and on the space coordinate z) a quasi-monochromatic 
beam is defined by 

El (i) = e iuJot A z x (t) , E z (t) = e lUot A z y (t) (7) 

where A. z , A z are stationary with stationary correlation and power spectrum 
inside (— u)i, uj{) . These properties remain whatever the chosen coordinates axes. 
It is reasonable to say that a wave is (purely) monochromatic when A. z x , A z verify 
a relation in the form (for some constant k and real B z (t)) 

A z x (t) = k cos B z (t) , A z y (t) = k sin B z (t) 

B z (t) represents the orientation and ke lu}at the complex amplitude of the field. 
This means that its (complex) amplitude is purely monochromatic in the usual 
sense. When B z is degenerate (B z (t) does not depend on t), the wave is po- 
larized (at z). In other cases, the wave is partially polarized. When E z and E 2 , 
have same power spectrum and are uncorrelated, the wave is unpolarized at z. 
Both properties are true in any orthogonal system. In optics, beams are often 
quasi-monochromatic and partially polarized but it is not always true, partic- 
ulary in astronomy and communications. Polarized and unpolarized beams are 
convenient idealizations. 

The concept of quasi-monochromatic wave is often bad-fitted even in the 
light domain. For instance Wolf-Rayet stars and B.L Lacertae have lines of 
relative width larger than few percents. Idem for LEDs (light-emitting diodes) 
with width larger than 10%. A more general model has to be taken in these 
cases. 

We have seen that the B.L law for one-component beams can be proved 
using a decomposition of media by LIF in series. The aim of this paper is to 
generalize the B.L law in the two-components case, using signal theory and 
modelling media as more general circuits. Inputs and outputs of these circuits 
are stationary processes which represent the components of the field. In the 
following section we consider that each component of the field at a distance z is 
the sum of two LIF outputs. The field components at the origin point (z = 0) 
are the inputs of these LIF. This model allows to determine the shape of the 
LIF characteristics, i.e the B.L law for bi-dimensional beams. 

Because two inputs and two outputs, 2x2 matrices of filters complex gains 
will be defined. In the years 1940, R. C. Jones had developed a "New Calculus 
for the Treatment of Optical Systems" 7 . It was based on 2x2 matrices which 
act in the time domain on purely monochromatic waves. Though formally re- 
sults of Jones are very close to formulas of this paper, they do not address the 
same objects. Actually Jones papers do not mention Beer neither Lambert or 
Bouguer, the pioneer of this topic. 

2 Two-dimensional Beer-Lambert law 

We deal with a beam which propagates in direction Oz of the orthogonal tri- 
hedron Oxyz. The electric field E z at time t is defined by its components 
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E* (t) , Ey (t) on axes Ox and Oy at distance z of the origin O. We assume that 
the medium between u and u + z is defined by a set of 2x2 "scattering matrice" 
H 2 independent of u 



I I Z l I Z 

TJZ TJZ 
11 Tin 



where the H^ k (to) depend on the frequency uj/2tt and are complex gains of LIF 
Hj k such that 



E%+* (t) = Hf 1 [E«] (t) + Hf 2 

e^ +z {t) = m x [e«] (t) + n z 22 



(t) 
(t) 



(8) 



Equivalently the electric field E z+ " at z + u is linearly dependent on its value 
E u at u [H] . The linearity is expressed by the set of the LIF Hj k which depend 
only on the medium. It is an obvious generalization of what was explained for 
the one-component waves propagation. Filters parameters depend only on the 
thickness of the considered medium. However results can be obtained without 
this hypothesis. Figure 1 shows equivalent circuits of (8) . 

Filters in series lead to multiplication of complex gains and filters in parallel 
to addition. The figure 2 summarizes the following equality 
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H*H< 
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(9) is equivalent to 
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We assume that the derivatives hP- 
(10) can be written as 



ok 



d_TT0 



H L are finite. The first equation of 



dz ±J -jk 



TJZ + U 



mi 



Hn - 1 



TJU 
TTZ H 21 



Obviously we have = H° 2 = 1 and H° 2 — — 0. When u 
(similar operation is done in the other equations) 



0, we obtain 



/if 2 = H\\h 12 

"21 — ■ n 2l' t ll 
^22 = ^21^12 " 



■ H[ 2 h 21 

- U{ 2 h 22 

- H 22 h 21 
H 22 h 22 . 
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The differential system can be split in two subsystems (equ.1+2 and equ.3+4). 
We assume that 



lim H 







(12) 



because any wave is evanescent in a passive medium. Two cases can be high- 
lighted following the (complex) eigenvalues Ai, A2 of the matrix 



"ll "-21 
^12 ^22 
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Solutions are in the form following that Ai ^ A 2 or Ai = A 2 = A 

1 jk 



or H jk = { c 3kiz + c jk2 )e Xz . 

The eigenvalues cannot cancel (or corresponding coefficients cancel), because 
(12) . Taking into account the initial conditions 

H° n = ff 2 ° 2 = 1 and H° 12 = H° 2l = (13) 



leads to two cases 
Case 1: Ai ^ A 2 

By identification with (11) we obtain 



n n 


= ae XlZ 


+ (1- 


- a)e X2Z 


H z 
n 12 


_ -fc? a 


(e XlZ 


- e X2Z ) 


a 21 


n 21 




- e X2Z ) 


A2 — Ai 


(e XlZ 


n 22 


= (1 - a 


)e XlZ 


+ ae X2Z 


a = 


A2-/1?! 

A2 — Ai 







(14) 



where the Aj,/i° fc can depend on frequency uj/2tt but are independent of z and 

X 1 = l(h° n +h 22 + ^pe i ^) 

A 2 = \ {h\ x + h° 22 - V/^ 2 ) (15) 
A=(h° 11 -h 2 > 2 ) 2 + 4h° 12 h° 21 =pe ie . 

The eigenvalues have real parts strictly negative (to fulfill the condition (12)). 
Case 2: Ai = A 2 

The solutions are given by the equalities 

H{ x = (az + 1) e Xz , H{ 2 = h° 12 e Xz 

H z 2 = {-az + 1) e Xz , ^ = h° 21 e Xz (16) 
a = 



"'11 ""ii \ ' t n~T"'22 

A — ^ . 



The case — h 22 =/= 0, h\ 2 — /i 21 — leads to the usual B.L law. These 
equalities are verified in any system of coordinates. Components evolve inde- 
pendently, with same attenuation and celerity. This corresponds to a medium 
with all possible properties of symmetry. 

In all cases, the real part of eigenvalues different from have to be negative 
for passive media which weaken waves. Moreover, the solutions are only matched 
to equations (11) , (12) , (13) . 

3 Examples 

In examples, we assume that the parameters h® k (w) are constant on spectral 
supports of inputs E°,E". 
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3.1 Example 1 

Let assume that 

£ h° 22 and h\ 2 = h° 21 = 0. 

We are in the case 1 with 

Hi x = e h °^ z ,H z 2 - e h °* z , H{ 2 = = 0. 

This means that a beam e lult polarized on Ox is transmitted with weakening 
exp[7\L z ] an d delay X [ftJJ z/oj. If polarized along Oy, the weakening is 
exp[ft [h° 22 ] z] and the delay is X [h° 2 ] z/cu. When TZ [h^] /TZ [h? 22 ] < 1, the 
first component disappears before the second component, independently of the 
values of X [h^] and X [h 22 ] which define the refraction indices of the medium. 
Then we can give to hi 1} h 22 values fitted to a dichroic material. More generally 
when 

££(*) = e h »*E%(t), E*(t) = e h °**E° y (t) 

both components evolve independently and the usual B.L law is verified for 
each component. The power P z at z becomes (we have assumed that the h® k 
are constant with respect to frequency) 



where at =E 



\E°At) 



\E° v (t) 



We retrieve the usual result when 



TZ [hii] = TZ [h 22 \ . In other cases usual B.L law is no longer valid because the 
two terms in P z have different behaviors. 

3.2 Example 2 

Let assume that we are in the case 1 (two distinct eigenvalues Ai, A2 different 
from 0) and that we deal with a monochromatic wave polarized along Ox: 

E° x {t) = e i »° t , E° y (t)=0. 

By (8) we have 

E^-WJJES], E*=?^[E°]. 
E° is transformed in E z defined by 




X2 — X1 



_ h 21 _ e X 2 z \ gi"0* 

A2— Ai V / 



where parameters can be complex. At z, both components are weakened and 
delayed through two terms functions of z (e XlZ and e X2Z ) and not one as in the 
usual B.L law. 

1) When h 21 = we have 

E% (t) = e h °^ z+tuJat , E z (t)=0. 

E z is polarized along Ox with TZ and X [ft-ii] as parameters of weakening 
and of delay. P z = e 2zK Ki] has the shape (5) of the usual B.L law. 
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2) Except when h 21 = 0, the monochromatic wave E z is no longer polarized. 
Assume that h° 21 ^ 0, h° 12 = (which implies ^ 0, h\ 2 ^ 0, ^ h% 2 ). From 
(14) 

J E z (t) = e h °i z+iuot 

I E y{ t ) = T^{ eh ° llZ - eh ° 22Z ) eiW0t - 

Even when 1Z = Tl [h 22 \ it is impossible to have P z like (5) . Consequently 
we understand using symmetries that usual BL law can be true for a particular 
polarization, and untrue for other polarizations. 

3.3 Example 3 

1) The case 

h° 12 = -h° 21 , h° 11 =hl 2 ^0 (17) 
is particularly interesting. We have 



\i = h\ 1 +ih\ 2 , \ 2 = h° 11 -ih a r 



12 



and (14) becomes 



H z 12 = -H z 21 = -i(e^ z -e^ z ). 



Now we assume that the electric field E° is monochromatic and polarized at 
angle <fi with respect to Ox 

El (t) = e iLJot cos (j) 
E° (t) = e l " ot sin 4>. 



From (8) and after elementary algebra 



E z (t) = e iuJot+h °i z cos (4> - zh° 12 ) 
E z (t) = e iwot+h °i z sin (0 - zh\ 2 ) . 



(18) 



For real hf 2 , the electric field E z at z is polarized at the angle (</> — zh\ 2 ) . Then 
a monochromatic wave at the frequency u!q/2it is rotated by the angle —zh\ 2 
and attenuated by exp {zR, [^ii])- Moreover —-^1 [^11] i s the time spent by 
the wave between O and z which shows that the medium refraction index n (uq) 
is equal to (c is the light velocity in vacuum) 

n(co ) = -—l[h° 11 \. 

The rotation angle is independent of the polarization angle <j). Because P z = 
e 2zTz[h 11 ] usua j b l i aw j s obeyed for the amplitude and the power, but not for 
components. 

2) When h\ 2 is not real (18) is developed in 

E x (*) = e*"°* +/l ii z (cosacosh6 + isinasinh6) 
E z (t) = e iuot+h ii z (sinacoshb - icosasinhb) (19) 



a = <j> - zR \h\ 2 \ , b = zl [h 
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Consequently ~E Z has two components polarized in orthogonal directions of an- 
gles a and (a + f ) with respect to Ox and of "height" Bf and B 2 such that 

Bf = e »[ fc ?i]* cosh (zl[h° 12 \) 
Bl = -^ e TC Kl] z sinh (zl [h° 12 ] ) (20) 
= tanh | zl [hi 2 ] \ ■ 



increases from to 1 when z increases from to oo. The main component 

is rotated by — zlZ \h\ 2 \ and attenuated by exp (zlZ [ftii]) cosh b. In the same 
time a second component appears which is orthogonal to the first component. 
The ratio of "heights" increases with z and tends towards 1. The power P z at 
z is given taking into account othogonality of components 

P z = e 2zTC Ki] cosh (2zl [h° 12 ]) 

which shows that the usual B.L law is not valid (except in the case of real h® 2 
which leads to P z = e 2zK [' 111 ]). The power of each component W x , W y does not 
follow the usual B.L law. For instance, for the power of the component 

pz = ^[fc?,] ( cog 2 ^ _ zU ^OJ) + sinh 2 ^0J)) _ 

The same remark is true for components in the directions a and (a + ^) with 
respect to Ox. 

3.4 Example 4 

We assume that E° is a quasi- monochromatic beam (see section 1) 
El (t) = e^Al it) , E° y (t) = e-ot^o {t) 

The h°- k (uo) are constant with respect of uj on the beam spectrum. If the con- 
ditions (17) are fulfilled, the components E°,EjJ are split by the medium in two 
orthogonal parts. With respect to Ox'y' defined by 

(Ox,Ox') = (Oy,Oy') = -zn [h° 12 ] 

the beam at z verifies using (19) and (20) 

f El, (t) = e^J' [A° x (t) cosh (zl [h° 12 ] ) + A y (t) sinh (zl [h\ 2 ] )} 
\ E* yl (t) = e R K^ [A y (t) cosh (zl [h° 12 ]) A® (t) sinh (zl [h° 12 ])] . 

We deduce the power P z defined by (6) 

\ P z = e 2K Ki] 2 [p Q cosh (2zl [h° 12 ] ) + sinh (2zl [h° 12 ])] 
\ e = 2l{E[A° x (t)A° y *(t)]}. 

Whatever the polarization state of the beam at z = 0, the usual B.L law is 
verified if and only if h\ 2 is real, i.e when the effect of the medium is a rotation 
(added to a weakening). Whatever h\ 2 , we have 9 = for instance when E° is 
polarized (E[..] is real) or unpolarized (E° and E° are uncorrelated) . 
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3.5 Example 5 

We assume that E° is a quasi-monochromatic beam like (7) and 

h° -h° h° - h° 

n ll ~ ;t 22i ""12 — "21 

i.c transfers between the coordinates are symmetric. From (8) we have 

#11 = #22 = eh ° lZ cosh ^ 2 z 
#12 = H 2i = e ' l ° lZ smhh° 12 z. 

However we remark that these relations are not maintained in other reference 
systems. For instance, in Ox'y' with (Ox, Ox') = <fi we have (the fey are the 
new parameters) 

fcOj = hl-L + /i? 2 sin2c?!> 
fc 22 = ^?i ~ ft-? 2 sin20 
k° 2 = fcOj = h° 2 cos 20. 

Elementary algebra leads to 

P z = e^Kh [p cosh (2zK [h° 12 ] ) + 8' sinh (2zft [/i° 2 ] )] 
0' = 2fc{E[4(*K*(*)]} 

which proves that the usual B.L law is not true apart from particular cases. 

4 Conclusion 

The Beer-Lambert law (actually due to P. Bouguer around 1729) was firstly 
used to measure the concentration of solutions . It addresses the problem of 
concentration measurement of some kind of molecules in a liquid. In equation 
1, we have 

a = k(uj)a (21) 

where a is the concentration and k (w) is a wavelength-dependent absorptivity 
coefficient, k is deduced from a measurement of the attenuation for a known 
value of a. The property of linearity with respect to the distance is due to Lam- 
bert and the linearity with respect to the concentration of absorbing species in 
the material was highlighted by Beer. The BL law intervenes in wave propaga- 
tion to explain together the attenuation and the dispersion whatever the crossed 
medium. 

A version of the Beer-Lambert law addresses the power as a function of the 
medium thickness, whatever the nature of the wave, acoustic or electromagnetic. 
The result has the form A (z) — A (0) e~ az where a is a function of the medium 
and of the frequency lo/2tt. Very often a is a power function of uj. Its estima- 
tion has numerous applications in medecine |13) . |16) . Also, chemistry uses the 
measurements of a because it is a linear function of the number of particles 
imbedded in the medium. The B.L law can be wrong for instance when multi- 
paths in fibres or in case of too strong transmitted powers [1] , [2] . Moreover the 
medium and the devices can be sensitive to polarization state. For instance the 
B.L law can be untrue in the case of birefringence for the medium and when 
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devices select only one component of the field. However a generalization is pos- 
sible studying separately both components of the electric field which defines an 
electromagnetic beam. 

The B.L law is easily proved modelling a medium thickness as a linear in- 
variant filter (LIF) where input and output show the evolution of the quantity 
of interest (for instance an amplitude or a power). In this paper we study the 
evolution of two quantities, the components of the electric field of an electromag- 
netic beam. To take into account interactions between components, a piece of 
medium is modelled by four LIF. We assume that the crossings of two successive 
medium pieces are independent events. This hypothesis suffices to determine 
the shape of the LIF complex gains. They are defined by four parameters h^, 
^12 i ^21 > ^-22 which depend on the medium and which may depend on the fre- 
quency. This model generalizes the Jones matrices used in the deterministic 
monochromatic beam to stationary random beams. As soon explained, Jones 
papers do not contain comments about the BL law. Examples of section 3 show 
that the set of parameters can be fitted to realistic situations. 
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